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Abstract 

What is the probability that a random triangle is acute? We explore this old question from a modern viewpoint, 
taking into account linear algebra, shape theory, numerical analysis, random matrix theory, the Hopf fibration, 
and much much more. One of the best distributions of random triangles takes all six vertex coordinates as 
independent standard Gaussians. Six can be reduced to four by translation of the center to (0, 0) or reformulation 
as a 2x2 random matrix problem. 

In this note, we develop shape theory in its historical context for a wide audience. We hope to encourage 
others to look again (and differently) at triangles. 

We provide a new constructive proof, using the geometry of parallelians, of a central result of shape theory: 
triangle shapes naturally fall on a hemisphere. We give several proofs of the key random result: that 
triangles are uniformly distributed when the normal distribution is transferred to the hemisphere. A new proof 
connects to the distribution of random condition numbers. Generalizing to higher dimensions, we obtain the 
“square root ellipticity statistic” of random matrix theory. 

Another proof connects the Hopf map to the SVD of 2 by 2 matrices. A new theorem describes three similar 
triangles hidden in the hemisphere. Many triangle properties are reformulated as matrix theorems, providing 
insight to both. This paper argues for a shift of viewpoint to the modern approaches of random matrix theory. 

As one example, we propose that the smallest singular value is an effective test for uniformity. New software is 
developed and applications are proposed. 
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1 Introduction 


Triangles live on a hemisphere and are linked to 2 by 2 matrices. The familiar triangle is seen in a different light. 
New understanding and new applications come from its connections to the modern developments of random matrix 
theory. You may never look at a triangle the same way again. 

We began with an idle question: Are most triangles acute or obtuse ? While looking for an answer, a note was 
passed in lecture. (We do not condone our behavior!) The note contained an integral over a region in M®. The 
evaluation of that integral gave us a number - the fraction of obtuse triangles. This paper will present several other 
ways to reach that number, but our real purpose is to provide a more complete picture of “triangle space.” 

Later we learned that Lewis Carroll (as Charles Dodgson) asked the same question in 1884. His answer for the 
probability of an obtuse triangle (by his rules) was 


8 - 
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TT 


0.64. 


Variations of interpretation lead to multiple answers (see Enn] and their references). Portnoy reports that in the 
first issue of The Educational Times (1886), Woolhouse reached 9/8 —4/7r^ « 0.72. In every case obtuse triangles are 
the winners - if our mental image of a typical triangle is acute, we are wrong. Probably a triangle taken randomly 
from a high school geometry book would indeed be acute. Humans generally think of acute triangles, indeed nearly 
equilateral triangles or right triangles, in our mental representations of a generic triangle. Carroll fell short of our 
favorite answer 3/4, which is more mysterious than it seems. There is no paradox, just different choices of probability 
measure. 

The most developed piece of the subject is humbly known as “Shape Theory.” It was the last interest of the first 
professor of mathematical statistics at Cambridge University, David Kendall inmni- We rediscovered on our own 
what the shape theorists knew, that triangles are naturally mapped onto points of a hemisphere. It was a thrill to 
discover both the result and the history of shape space. 

We will add a purely geometrical derivation of the picture of triangle space, delve into the linear algebra point 
of view, and connect triangles to computational mathematics issues including condition number analysis, Kahan’s 
[315] accuracy of areas of needle shaped triangles, and random matrix theory. 

We hope to rejuvenate the study of shape theory! 



14 PILLOW-PROBLEMS. 


57. (aSi So) 

In a given Triangle describe tKree Squares, whose bases 
shall lie along the sides of the Triangle, and whose upper 
edges shall form a Triangle; 

(i) geometrically; (2) trigonometrically. [27/1/91 

58. (^ 5 , 83 ) 

Three Points are taken at random on an infinite Plane. 
Find the chance of their being the vertices of an obtuse- 
angled Triangle. [20/1/84 


Figure I: Lewis Carroll’s Pillow Problem 58 (January 20, 1884). 25 and 83 are page numbers for his answer and his 
method of solution. He specifies the longest side AB and assumes that C falls uniformly in the region where AC 
and BC are not longer than AB. 
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1.1 


Random Angle Space 


( 0 , 180 , 0 ) 



Figure 2: Uniform angle distribution: angle 1+ angle 2+ angle 3 = 180°. 

The simplest model of a random triangle works with random angles. We mention it as a contrast to the Gaussian 
model that is more central to this paper. Figure shows a picture of angle space in K^. It is a “barycentric” picture 
in the plane for which angle 1 + angle 2 + angle 3 = 180°. A random point chosen from this space has a natural 
distribution, the uniform distribution, on the three angles. We then reach a simple fraction: 3/4 of the triangles are 
obtuse. 

The normal distribution on vertices also gives the fraction 3/4. This result is much less obvious. We are not 
aware of an argument that links the “angle picture” with the normal distribution, though it is hard to imagine that 
anything in mathematics is a coincidence. Nonetheless, the normal distribution on vertices gives a very nonuniform 
distribution on angles. We explore this further in Section 3.6. 

In the next section we will discuss the “natural” shape picture of triangle space that is equivalent to the normal 
distribution. It seems worth repeating a key message: the uniform angle distribution is not the same triangle measure 
as the Gaussian distribution yet obtuse triangles have the same probability, 3/4. 

1.2 Table of Random Triangles 

Before computers were handy, random number tables were widely available for applications. On first sight, a booklet 
of random numbers seemed an odd use of paper and ink, but these tables were highly useful. In the same spirit, 
we publish in Figure a table of 1000 random triangles. Only the shapes matter, not the scaling. You might try 
to count how many are acute. The vertices have six independent standard normally distributed coordinates. Each 
triangle is recentered and rescaled, but not rotated. 
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Figure 3: 1000 Random Triangles (Gaussian Distribution): Most triangles are obtuse. 

1.3 A Fortunate “Optical Illusion” 

This subsection moves from 1,000 to 50,000 random triangles. The assumptions are the same: cci,j/i, 0 : 2 , 2 / 2 , 2 : 3 ,j /3 
are six independent random numbers drawn from the standard normal distribution (mean 0, variance 1). These six 
numbers generate a random triangle with vertices {xi,yi), (x 2 , 2 / 2 ) and (xs^ys). From the vertices, we can compute 
the three side lengths, a, 6 , c. And since we do not care about scaling, we may normalize so that a^ + 6 ^ + c^ = 1. 

Instead of drawing the actual triangle, we represent it by the point (a^, P, P) in the plane x + y + z = 1. The 
result for many random triangles appears in Figure]^ 

The first curiosity is that the collection of points forms a disk. This is the triangle inequality, though the 
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connection is not obvious. Points on the outer circumference represent degenerate triangles, with area 0. 

A second curiosity concerns right triangles. The points that represent right triangles land on the white hgure, an 
equilateral triangle inscribed in the disk. 

A third curiosity, not visible in the picture, and also not obvious, is that one quarter of the points land in the 
acute region, inside the white equilateral triangle. Each of the three disk segments representing obtuse triangles also 
contains one quarter of the points. 

A fourth curiosity, perhaps the most important of all, is the particular density of points (triangles) towards the 
perimeter of the disk. It is not difScult to imagine, in the spirit of many familiar optical illusions, that one is looking 
straight down towards the top of a hemisphere. The three white line segments are semicircles on the hemisphere, 
viewed “head on”. 



Figure 4: Points represent triangle shapes, where (a^, 5^, (?) is drawn in the plane x + y + z = 1. Right triangles fall 
on three white lines like = 1/2. 

Indeed, the points are uniformly distributed on a hemisphere. This fact was known to David Kendall, and is a 
major underpinning of the subject of “shape theory.” We discovered it ourselves through this picture. Each area is 
1 /4 of the area of the hemisphere, as Archimedes knew. 

Grade school geometry emphasizes “side-side-side” as enough to represent any triangle, so what does height on 
the hemisphere represent ? The answer is simple: height represents the area. The Equator has height zero, for 
degenerate triangles. The North Pole, an equilateral triangle, has maximum height. Latitudes represent triangles of 
equal area. 

There is more. Triangles with one angle specified form small circles on the hemisphere going through two vertices 
of the white figure. Indeed it is good advice to take any triangle property and consider what this hemisphere 
representation has to say. 
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1.4 Two by Two Matrices: Turning Geometry into Linear Algebra 

A triangle may be represented as a 2 by 3 matrix T whose columns give the x,y coordinates of the vertices. One 
way to remove two degrees of freedom is to translate the first vertex, say, to the origin. A more symmetric way to 
remove two degrees of freedom is to translate the centroid to the origin. There is something to be said for treating 
all vertices and also all edges symmetrically. The classical law of sines does this, as we will do, but not the law of 
cosines. Hero’s formula for area is also symmetric but ^ base x height requires the choice of a base and hence is 
not symmetric. 

Let A be the 2x3 matrix of a reference equilateral triangle centered at the origin. This means that each of the 
three columns of A has the same euclidean length. A 2 x 2 matrix M transforms the triangle by taking the vertices 
of the equilateral triangle to the vertices of another triangle centered at the origin. As illustrated in Figure]^ below, 
every 2x2 matrix M may be associated with a triangle with zero centroid through the equation T = M/S.. If the 
matrix M is random, then the associated triangle T is random. 



Figure 5: Equivalence between 2x2 matrices M and triangles. 


An origin centered triangle may alternatively be represented by its edge vectors which we may place in a 2 x 3 
matrix E whose columns add to the zero column (0,0)^. We may again let A be the notation we use, this time 
for the edge vector matrix for an equilateral triangle. Whether A is thought of as the vertices of an origin centered 
equilateral triangle or the edges, the important property is that the columns of A are all of equal length. The vertices 
are equidistant from (0,0), the edges have equal lengths. 

If M is a random 2x2 matrix, then E = MA produces a random triangle whose columns are triangle edges. 
The columns of A sum to (0,0)^, hence the columns of E sum to (0,0)^, so E encodes a proper closed triangle. The 
edge lengths are the square roots of the diagonal elements of E^E. The vertex view would have produced T, whose 
three column vectors originate at the centroid, giving less access to the edge lengths. 

Section 2.2.1 takes a close look at the choice of the Helmert matrix A = ^ 1/^ 1/^ 2/i/6 ) 

reference. 

Consider as an example, the “45,45,90” zero centered right triangle with vertices (—2, —1),(1, —1),(1, 2) and edge 
lengths 3-\/2,3,3. We form the matrix T = ^ ^ 1 2 ) ' Helmert matrix as A, the corresponding 

My = TA^ for the vertex view is My = — ^ translation of the triangle produces the same 

My. 

An edge view requires differences of columns of T. We may transform between 0 centered vertices and edges with 
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the equations 



E = T 


These matrices are pseudoinverses as the (1,1,1) direction is irrelevant. 

/ _3 3 Q 

The 3-\/2, 3,3 triangle is then represented as E = 


VTs 0 


-3 0 3 


. The corresponding edge view matrix is Me = 


EA^ = - 

It is possible to check that My = M, 


1/2 V3/6 

-v^/6 1/2 

This paper will take the edge view unless noted otherwise. 


always holds. 


2 Triangular Shapes = Points on the Hemisphere 

Now we come to the heart of the paper. Every triangular shape with ordered vertices is naturally identified with 
a point on the hemisphere. “Random triangles are uniformly distributed on that hemisphere.” We know of three 
constructions that exhibit the identification: 

• Complex Numbers 

• Linear Algebra: The Singular Value Decomposition 

• High School Geometry 

The linear algebra approach, through the SVD, generalizes most readily to higher dimensional shape theory. It also 
invites applications using Random Matrix Theory which has advanced considerably since the early inception of 
shape theory. (See Sections 3.2 and 4.1 of this paper for connections to Random Matrix Theory including random 
condition numbers and uniformity tests.) Yet another benefit is the connection to the Hopf fibration. 

The geometric construction has two benefits not enjoyed by the others: 

1. It can be understood with only knowledge of ordinary Euclidean geometry 

2. It reveals the triangular shape and the point on the hemisphere in the same picture. 

2.1 Complex Numbers 

A generic triangle can be scaled, rotated, and reflected so as to have vertices with complex coordinates 0,1, (/• There 
are six possibilities for ( in the upper half plane, corresponding to the six permutations of the vertices. One can 
consider C in the lower half plane, to pick up six more possibilities. 

The usual correspondence with the hemisphere is that the stereographic projection maps (/ from the upper half 
plane to the hemisphere mm- 



2.2 Linear Algebra 


2.2.1 The Helmert matrix and regular tetrahedra in n dimensions 


The Helmert matrix A appears in statistics, group theory, and geodesy. It is a particular construction with orthonor¬ 
mal rows perpendicular to (1,1,1): 


f 1/^2 -1/^2 0 \ 

1 / V 6 1/^6 - 2/^6 J ■ 


The columns contain the vertices of an equilateral triangle. They are also the edges of a similar equilateral triangle. 
The operator view is that is an isometry (distances and angles are preserved) from to the plane x + y + z = 0. 
This means that AA^ = I 2 and A^A = — J 3 / 3 , where J 3 =ones(3) is the all-ones matrix. 

A useful property related to the duality between vertices and edges is that 


A 


1 

-1 


-1 

1 



= 73 ^ 


- sin f 

cos f 
6 


A. 


The left hand side takes the vector differences of the columns of A (the vertices of the equilateral triangle). These 
differences are the edges. We may place them as vectors at a common origin. The right hand side is -s/S times a 
rotation matrix times A. 

The -s/S represents the familiar fact that the edge lengths of an equilateral triangle are -s/S times the distances of 
the vertices from the centroid. Also if you rotate an edge counterclockwise by 7r/6, you get the direction of a vertex. 

The generalization of Helmert’s construction to K" is the n — 1 x n matrix A„. It has orthonormal rows, and the 
zeros make it lower Hessenberg: 


/ 



/ 1 -1 0 . 

1 1 -2 . 

• ® \ 
. 0 

V 

yj{n-l) -n ] 


Iv 1 11. 

■ -{n-1) j 


A transposed (and negated) form of that last matrix is available in the statistics language R as contr. helmert (n). 
The name indicates a matrix of “contrasts” as used in statistics. 

The square (and orthogonal) Helmert matrix adds a first row with entries (1,... ,l)/y^n. It is obtained in 
MATLAB as an orthogonal test matrix with the command gallery('orthog’ ,n,4). 

The rectangular A„ contains the vertices of a regular simplex from the “column view.” When n = 4 this is the 
regular tetrahedron with four vertices in Our equilateral triangle A is A3. (The “row view” consists of n — 1 
orthogonal rows each perpendicular to (1,1,..., 1) in K".) Again A„A^ = /„_i and A^A„ = /„ — J„/n. 


2.2.2 The SVD : How 2x2 matrices connect triangles to the hemisphere 

Let M = U'EV"’" be the Singular Value Decomposition of a 2 x 2 matrix M. As shown in Figure with the edge 
viewpoint, the columns of A = MA are ordered edges of a triangle. We take a few steps to make the shape and the 
SVD unique: 

• Scaling: We may assume that the squared entries of M have sum 1. The diagonal matrix E then has 1 = 
ctJ -I- (t|. The triangle edges then have sum of squares 1 since tv{E'^E) = tr(A^M^MA) = tr(M^MAA^) = 
tr(A/^M) = 1. 
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• The orthogonal factor U in the SVD is unimportant, since M and U correspond to the same triangle just 
rotated or reflected. 

I COS 0 _ sin 0 \ 

• To make the SVD unique, assume that ai > a 2 > 0 and that V = [ . „ „ with 0 < 9 < tt. There 

~ ~ sm0 costl J ~ 

is a singularity in the SVD when M = ^12- In this case V is arbitrary. 

We associate V and E with a point on the hemisphere of radius ^ : The longitude is 26 and the height is CTitT 2 - The 
latitude is thus asin(2cricr2). From 1 = af + we have 0 < CTicr 2 < The singularity of the North Pole (every 
angle is a longitude) is consistent with the singularity of the SVD at :^^ 2 - 

The height (71(72 is also (k + where k = <7iI<J2 > 1 is the condition number of M. The best conditioned 

matrix (k = 1) is at the North Pole, and corresponds to the equilateral triangle. The ill-conditioned matrices 

with K = cx), and at height 0, are on the equator. They correspond to degenerate triangles, with collinear vertices. 
We see our first hint of the link between “ill-conditioned” triangles and “ill-conditioned” matrices. 

One can now go either way from triangles to the hemisphere through the 2x2 matrix M: 

• Hemisphere to matrix to triangles: Start with a point on the hemisphere. Create M = 'EV'^ from the latitude 
asin(2(Ti(T2) and longitude 29. The edges of the triangle are the columns of MA. 

• Triangles to matrix to hemisphere: Start with a triangle centered at 0 whose squared edges sum to 1. Then 
M = (2 X 3 matrix of triangle edges)A^. The point on the hemisphere comes from the SVD of M by ignoring 
[7, taking the height from detM = (Ti(T 2 , and the longitude from 9. 

The matrix M is related to the preshape in [9]. In that view, any triangle can be moved into a standard position 
through A^. Our recommended view as shown in Figurej^is that the preshape is an operator mapping the equilateral 
triangle into a particular triangle. This viewpoint emphasizes the linear operator that transforms triangles. 


2.3 Formulas at a glance 


As a summary, one may sequentially follow these steps to derive the key formulas: 

Reference Triangle: edges from the three columns of A = 

Random Triangle: edges from the columns of E = MA. (Note M = EA"^ since AA^ = I.) 


f 1/^2 -1/^2 0 \ 

'v 1/V6 1/V6 -2/V6 J ■ 


(1) SVD ((71,(72,0) of Matrix M (||M|| = 1) : 

\ (72 J \ — sm.9 cos 9 J 

1 > (7i > (72 > 0, (7j -I- (7| = 1, 0 < 0 < TT, (t/ not needed.) 


(2) Triangle edges a, b, c: 

j 

= 1 , 6 ^ 


= diag(A’^M^MA), a-l-6>c, 6-|-c>a, c-l-a>6. 


(3) Hemisphere of radius 1/2 (coordinates A, (/ denote I(cos Acos(/>,cosAsin^, sinA): 

Latitude = A = asin(2(7i(72), Longitude = cj) = 29, Height = ^ sin(A). 0 < A < 7r/2, 0 < (j) < 2 -k. 


10 


(4) Disk of radius 1/2 (projection of hemisphere to polar coordinates r,(l) ): 

r = i cos(A). The angle ^ = 20 is the longitude of the point on the sphere. 0<r<l/2, 0<(/)< 27r. 

The area formulas for these descriptions: 

K = Area = = i^1 - 2{a^ + h^ + c^) = ^ = {k + 

Here are conversion formulas. The (2)—>-(4) and (4)—>-(2) blocks are developed in Section 2.3.1. 


From\To 

(1) SVD(M) 

(2) Triangle 

(3) Hemisphere 

(4) Disk 

(1) SVD(M) 

CTl, 0-2,0 


diag((MA)^(A/A)) : 
a2 = i(l — (cr| — cr|) cos20-|_) 
52 = 1(1 — {a\ — cr|) cos 20_) 
c2 = 1(1 — (cr| — (t|) cos 20) 

A =asin(2cricr2) 

</=20 

Height=AT-\/l2 

= aia2={K+K~^)~^ 

rsin(/)= {M'^M)i 2 
rcos0 = 

^ = \J\~ ^1^2 

= \{crl-al) 

^ = 20 

(2) Triangle 
a^, 52^ ^2 

1) Triangle—>■ disk 

2) then use (4) 1 


j_ /cosAsin(/A _ A ( “2 1 
Vg y cos a cos 0 J 1 2 ) 


(3) Hemisphere 
X,(j) 

cji = cos(A/2) 
CJ 2 = sin(A/2) 

1) r = cos(A)/2 

2) then use 1 


r = cos(A)/2 

(4) Disk 
r,(j) 

cr| = 1/2 + r 
an = ^li-r 

0 = <//2 

a2 = (1 — 2rcos(/+)/3 

52 = (1 — 2rcos(/_)/3 
c2 = (1 — 2rcos(/)/3 

A = acos(2r) 

4> = 4> 



0± = 0 ± 7r/3, K = cri/(T 2 , (/± = (/ ± 27r/3. 


Also useful is the direct mapping of M to the hemisphere in Cartesian coordinates: 


M = 


Mil Mi2 
M 21 M 22 


/ {Mh + Mil) - + MI2) 

2 (MiiMi 2 + M21M22) 
y 2IM11M22 — M21M12I 


/ (M^M)n - (M^M)22 
2{M^M)i2 
y 2|detM| 


^ cos A cos 
cos A sin (j) 
y sin A 


This formula may be derived directly or through the Hopf Map described in Section 2.5.5. 

A notebook in the programming language Julia is available which takes six representations and performs all 
possible conversions. The notebook provides a correctness check by performing all possible roundtrips between 
representations. (See http://www-math.mit.edu/~edelman/Edelman/publications.htm,) 


2.3.1 Disk r,cj) to triangle a^,b^,c^ directly 
An immediate consequence of 


a 

^2 


= diag((MA)^(MA)) = diag(A^HE2H^A) 
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= ^I + r 


is that 


This contains the quadratic form V 
One can readily derive 


1 1 +r I 


evaluated at the three columns of A. 


We prefer the geometric realization that A contains three columns oriented at angles tt/G, Gtt/G, Qtt/G with lengths 

■^2/3. The matrix ^ ^ ^ ^ rotates a vector at angle a clockwise by d, reflects across the x-axis, and then 

rotates counterclockwise by 0. The quadratic form thus takes a dot product of a vector at angle a — 9 with its 
reflection, i.e., two vectors at angle 2{a — 9), yielding cos(2a —2(/)) = — cos{2(j) — 2a + tt). Plugging in the three values 
7r/6, Stt/G, and Qtt/G gives the result with little algebra! 

As (j) runs from 0 to 27r, the points (sine/),cos(/)) trace out a circle in the plane. The transformation to the plane 
x + y + z = 0 is 


= — COS 0- 


This yields an alternative formula to transform the disk to the squared sidelengths: 


1 l+^/grA^ 


2.3.2 Triangle to disk r,(j) directly: a “barycentric” interpretation 

Apply A to the equation just above: 


A 62 


Then r and </> are the polar coordinates of 


A 62 = 


= vr 


where A = 


1/2 1/2 -1 
V3/2 -V3/2 0 


Barycentric Interpretation: To obtain the point on the disk, take an equilateral triangle of sidelength 1 (the 
columns of A) and find the point with barycentric coordinates {a^, 62, c2). 
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The “broken stick” problem: asks how frequently a, b, c can be edge lengths of a triangle, when their sum 

is the length of the stick. How often do a,b,c satisfy the triangle inequality? This has an easy answer 1/4 (again!) 
and a long history. 

We could modify the question to fix = 1. The answer now consists of those (a, b, c) inside the disk 

in Figure 7(a) as a percentage of the big triangle. The solution to our “renormalized” problem is that the fraction 
7r/-\/27 « 0.60 of all triples {a,b,c) can be edge lengths of a triangle. 

2.3.3 Triangles with fixed area 

In the four representations the following “sets of fixed area” are equivalent: 

• Matrices M with (||M|||. = = 1 and) determinant ^/VIK. 

• Triangles with (a^ + &^ + = 1 and) constant area K = ^(Hero’s formula). 

• The latitude at height \/l2K on the hemisphere. 

• A circle centered at the origin with radius r = — 48AT^ in the disk representation. 

To sweep through all triangles with area K, take r = ^\/l — 48K^ and let (f) = 29 run from 0 to 27r/3. This 
avoids the “rotated” triangles that cycle through a, b, c. 

Here are special triangles with area K (including approximations for very small K): 


• Right Triangles: 


K 

r 

(j) (or (/)±) 

squared sides 

K < 1/8 

r > 1/4 

±acos(-^) 

i and 


For area K < ^ (corresponding to r > |) there are six congruent right triangles. 

For K small, the right triangles have squared side lengths 8K^ + 128K^ + O(K^), and ^ — 128iF^ + 

0{K^). We can check that ^^/{8K^){l/2) « K. 


• Isosceles Triangles (third side smaller): 


K 

r 

(j) (or (/±) 

squared sides = 1 

K < l/^/^ 

r < 1/2 

0 

Two equal sides:(l + r)/3 Third side:(l — 2r)/3 


For each K, three congruent isosceles triangles might be called the most acute of the acute triangles (furthest 
from the three white lines that represent right triangles.) As r —>■ i, this isosceles triangle approaches a right 
triangle. (Many physics computations use this fact.) The squared side lengths for K small are two of size 
i - 4K^ - 48K^ + 0{K^) and a tiny side 8K‘^ + 96 + 0{K^). 


• Isosceles Triangles (third side larger): 


K 

r 

(j) (or (/±) 

squared sides + &^ + = 1 

K < 1/V48 

r < 1/2 

7r/3 

Two equal:(1 — r)/3 Third:(1 + 2r)/3 


For each K, three congruent isosceles triangles might be called the most obtuse of the triangles (although they 
may be acute). The squared side lengths for K small are two of size | + 4K^ + 487F^ + 0{K^) and a third side 
l-8K^ + 96K^ + 0{K^). 
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• Singular Triangles: 


K 

r 

(j) (or (j)±) 

actual sides a,b,c 

0 

1/2 

any 

\/i 1 sin i(/I, y 11 sin i(/+1, y 11 sin i(/_ 1 


The longer side is the sum of the two shorter sides. 


• Nearly Singular Triangles: 

For K tiny, r = i - - lUK'^ + 0{K^). 

We conclude with area/angle formulas that will be useful in Section 2.5.3. They may have interest in their own 
right as they do not seem to appear in the famous Baker listings of area formulas from 1885. [T]. 

Lemma 1. The following formulas relate side lengths a,b,c, area K, and angles A,B,C: 

tan A = A:K/{tP‘ + — a^) 

tan B = 4K/(a^ + — b^) 

tan C = 4K/(af + b^ ~ c^)- 

We can derive tan A from the law of cosines: cos A = —a^)/26c and the well known formula sin A = 2K/bc. 

We apply Lemma 1 with = 1 in the form 

tan A = 4iF/(1 — 2a?) 
tan B = 4K/\l - 2b^) 
tan C = 4iF/(l — 2c^). 

If angle A is held fixed, there is a linear relationship between and K. This means that a circular arc on the 
hemisphere represents all triangles with angle A. This will be further discussed in Section 2.5.3. Three special arcs 
(the three white lines in Figure 3) represent all right triangles with A or B or C equal to 7r/2. 

2.4 The Triangle Inequality and the Disk Boundary 

The razor sharp triangle inequality shows up as the smooth disk in the plane x -\- y -\- z = 1. We saw this in Section 
1.2 experimentally. Careful readers may have thought about this disk throughout Section 2.3, and some points are 
worth attention: 

• Area Viewpoint: The triangle inequality is equivalent to non-negative area. Hero’s formula 

K = i a/I — 2{a'^ + b‘^ + (A) implies -I- 6^ -I- < 1/2. This is the intersection of the sphere of radius l/\/2 

centered at the origin with the plane x-\-y-\-z = 1 whose distance from the origin measured from (1/3,1/3,1/3) 
is l/-\/3. The radius of the circle of intersection is then ^1/2 — 1/3 = \/l/6. Equation ([^ takes us to the disk 
with r = 1/2 through the factor of -^3/2. 

For general a, b, c (no restriction on -I- 6^ -|- c^). Hero’s formula is 

16A^ = (a -I- 6 -I- c)(—a + b + c)(a — b + c)(a + b — c) = (a^ + b^ + — 2(0"^ -I- 6"^ -I- o'*). 

This formula has three factors whose positivity amounts to the three triangle inequalities. The inequality 
(x + y + z)^ — 2(x^ + y^ + z^) > 0 produces a cone through the origin which intersects the plane x + y + z = 1 
at the aforementioned circle of radius l/-\/6- 

• Matrix Viewpoint: The boundary of the triangle inequality is equivalent to | detM| = aia 2 = 0. 
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• Hemisphere Viewpoint: det(M) = 0 corresponds to the equator of the hemisphere (A = 0). 

• Disk Viewpoint: The conversion table yields r = cos(A)/2 < 1/2, where A = 0 corresponds to the disk 
boundary at r = 1/2. 

• Triangle Viewpoint: The sides of the singular triangles (with = 1) are expressed in the table 

before Lemma 1. 

2.5 Geometric Construction 

Here we exhibit what is marvelous about mathematics: The triangle is not only represented as a point on the 
hemisphere, but it can be constructed geometrically inside that hemisphere. The full construction will be described 
in a moment. Informally, the point on the hemisphere will be one vertex, and the base is inscribed in an equilateral 
triangle on the equatorial plane. 

2.5.1 Parallelians 

We must define the three parallelians of an equilateral triangle E, meeting at an arbitrary point P. 

Definition 2. (Parallelians) Through any point P, draw lines parallel to the sides of E (Figure^. The parallelians 
are the segments of these three lines with endpoints on the sides (possibly extended) of E. 

When P is interior to E, the three line segments are interior as well. They intersect at P. Otherwise their 
extensions intersect at P. 




Figure 6: Parallelians when P is interior or exterior to E. 

2.5.2 Barycentric coordinates revisited 

This section picks up where Section 2.3.2 left off, with a barycentric view of the point in the disk. We introduce big 
and little equilateral triangles whose vertices are the columns of 

A _ { V3/2 -V5/2 0 \ ^ . 

big 1/2 1/2 -1 ) 2 
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^little 


/ 0 

1/2 

1/2 

:=Abig 1/2 

0 

1/2 

V 1/2 

1/2 

0 

or equivalently 

^big 

= A] 


-v^/4 v^/4 0 

-1/4 -1/4 1/2 


1 . 

-2^big 



The columns of all have norm 1. The midpoints (columns of all have norm 1/2. 

We have seen in Section 2.4 that if the sides have = 1, then 



^little 


— 

— 6^ + (? 
+ 6^ — (? 


1 - 2a^ 

= %ttle| 1-2&; 

1 - 2c2 


falls inside a disk of radius 1/2. (Note when comparing Section 2.3.2, for convenience, we have reflected the x and y 
coordinates in A. This serves the purpose of allowing the little triangle to have a horizontal base.) 

In summary, the barycentric coordinates (a^,6^,c^) for the “big” equilateral triangle correspond to the (possibly 
nonpositive) barycentric coordinates (1 — 2a^, 1 — 2&^, 1 — 2c^) on the inverted “little” triangle. (Figure 7(a)) 

Figure 7 also includes the segments into which parallelians are sliced. In an equilateral triangle, the barycentric 
coordinates are in the same ratio as the distances to the sides in (d), which in turn are in the same ratio as the 
parallelian segments in (e). We note that the circle has radius 1/2, the little triangle has edges •\/3/2 and the big 
triangle has edges -v/S. In Figure 7(e), the lengths labeled a,/3,7 are exactly ^(1 — 2a^), ^(1 — 25^), ^(1 — 2c^). 


(a) (b) (c) (d) (e) 



Figure 7: (a) P = (a^,&^,c^) in the A^jg barycentric system, (b) P = (a,/3,7) in the A^j^^jg system. Also in the 

same ratio a : (3 : ^ are the areas (c), the perpendicular distances (d), and the parallelian segments (e). Explicitly 
a : /3 : 7 = (1 - 2a^) : (1 - 2&2) : (1 - 2c^). 

Endpoints for Parallelians through P in Barycentric Coordinates 
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<1 

^little 

p 

(a2,62^c2) 

(1 - 2a2,l - 262,1 - 2 c2) 

Parallelian 1 

(a2, 1/2 ,l/2-a2) 

(a2, 1/2-a2, 1/2) 

(1 - 2a2, 0 , 2a2) 

(l-2a2, 2a2, 0) 

Parallelian 2 

(1/2 ,62,1/2-62) 
(1/2 - 62, 62, 1/2) 

(0 ,1 - 262,262) 

(262, 1 - 262, 0) 

Parallelian 3 

(1/2 ,1/2-c2,c2) 

(1/2 - c2, 1/2 , c2) 

(0 , 2c2, 1 - 2c2) 

(2c2, 0,1- 2c2) 


2.5.3 Hemisphere to Triangle: Direct Construction 

Let S' be a point on the hemisphere corresponding to edge lengths a, b, c. A straightforward but unusual construction 
yields a triangle within that hemisphere whose sides are proportional to a, 5, c. In the construction below, we let E 
denote the equilateral triangle inscribed in the hemisphere’s equator as in Figure 7(a). The vertices of E in the xy 
plane are the columns of 

Geometry Construction Theorem: Let S be any point on a hemisphere of radius 1/2. Let P be its projection 
onto the base. The line through P parallel to the x-axis, intersects A at X and Y. The triangle SXY in Figure]^ has 
side lengths proportional to a, 6, c. 

Proof: Let uj = We know from the last sentence in Section 2.5.2 that XP=cj(f — h^) and YP=cj(f — a?'). 
Thus the parallelian length is XY=a;c^. The height SP is \/\2K = 2ujK from the table in Section 2.3. Thus the 
tangent of angle SXP is 4AT/ (1 — 26^) and similarly the tangent of angle SYP is AK/{I — 2a^). These agree with the 
tangents in Lemma 1, so we have constructed the triangle with sidelengths proportional to a, 6, c. 

The lengths are SY=ujac, SX=ujbc and XY=a;c^. This triangle inside the hemisphere is one of the three similar 
triangles that we now construct using parallelians. 

The collinear points P, X, Y in the xy plane are the columns of 


[P X Y]= A 


big 


62 


1/2 

1/2-c2 


1 / 2 -< 
1/2 


From a,b^c we computed P,X,Y in the equatorial plane in Figure 9 from this formula. We construct the 


hemisphere point S by changing P’s z-coordinate to 


i-\\pv- 
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xy plane 



Parallel to xz plane 


Figure 8: Top View (Left) and Front View (Right) 

Triangle SXY has sides proportional to a, b, c. 


2.5.4 The three similar triangles theorem 

Theorem 3 is interesting for several reasons. It seems unusual that in traditional geometry, one would construct 
triangles with vertex on a hemisphere, and opposite edges inscribed in an equilateral triangle that itself is inscribed 
in the hemisphere’s base. It also feels unusual that the three triangles would share a common line segment SP which 
represents three different altitudes in the shape. 

Theorem 3. S is on the hemisphere corresponding to a,b,c and P is its projection onto the base. Three triangles 
share SP as an altitude with vertex S and opposite edge one of the parallelians through P as in Figure 7(e). The 
three triangles thus created are similar with edges in the proportion a : b : c. 

Proof: The construction in 2.4.3, by symmetry, yields three triangles. The scaling factors 2uja and 2ojb and 2wc 
multiply a,b,c. The three similar triangles are illustrated in Figure 9. 
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Illustration of the hidden geometry inside triangle shape space revealed in Theorem 3: 
Figure 9: The three similar triangles theorem. 

Above: Skeleton: The three triangles sharing the altitude SP are similar! 

Below: 3d model: All triangles are similar; congruent triangles are color coded. 
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2.5.5 Matrix to Hemisphere: The Hopf Map 

Hopf discovered in 1931 (see [12]) that the hypersphere maps naturally to the sphere 5^. We think of points on 
as the 2x2 matrices normalized by ||M|||, = ^ Mfj = 1. The Hopf map comes directly from the Singular Value 
Decomposition M = 


Hopf(M) ■.M = U 


cos(A/2) 0 
0 sin(A/2) 


cos( 0 / 2 ) — sin(^/ 2 ) 

sin((/)/ 2 ) cos(^/ 2 ) 


/ cos A cos( 
cos A sin Cf 
^ sin A 


52 . 


Here the latitude is A G [— 7 r/ 2 , 7 r/ 2 ], with sign matching det(M). The |A| is determined by the singular values of 
M, while the longitude (j) comes from the right singular vectors of M. The rotation matrix U of left singular vectors 
gives the fiber of that is mapped to a single point on S'^. 

Remark: We have not seen the SVD in the usual definitions of the Hopf fibration. It may be the the most 
immediate way to get to the Hopf map. 

An explicit formula is 


Hopf(M) : 


Mn 

M 21 


Mi 2 

M 22 


[Ml, + Mil) - (MI2 + MI2) 
2 (MiiMi 2 + M21M22) 
2 {Mi\M 22 — M21M12) 


Our construction of triangle space, mapping M to the upper hemisphere, uses ^Hopf(M). The coordinate, 
z = I sin A, is always chosen positive. 


2.5.6 Rotations, Quaternions, and Hopf 

The bad news is that there is an ugly formula for the general 3x3 rotation matrix: 


+13^ 

Qs = I -2(0(5 - j3"f) 

2 ( 0:7 + 136) 


2 {ad + f3"f) 

q ,2 — /JS _|_ 72 — ^ 

—2(0/3 — 7(5) 


2(/3(5 — 07 ) 
2(o/3 + jS) 

o2 — /32 — 72 _|_ J2 


where 1 = o^ + /3^ + 7 ^ + This is a rotation about axis (/ 3 , 7 , 6 ) with rotation angle 2 acos(o). 

The good news is that given any Q 3 we can obtain o, /3, 7 ,6 through an eigendecomposition. We can then associate 
a 4x4 rotation matrix 


Q 4 — 


The key relationship that ties Q 3 to Q 4 will be crucial for Theorem 6 : 

Hopi^QiM) = (53Hopf(M). 

Q 4 M denotes the operation of flattening M to the column vector [Mu, M 21 , M 12 , M 22 ]^, applying ( 54 , and reshaping 
back into a two by two matrix. 


/ a 

-P 

5 

-7 \ 

P 

a 

7 

(5 

-5 

-7 

a 

P 

\ 7 

-5 

-P 

a / 
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Readers who prefer not to think about quaternions can safely skip ahead to Section 2.6. Hopf(Q4M) = 
(53Hopf(M), could be checked directly, for example by Mathematica. But quaternions are the better way to under¬ 
stand why the relationship holds. It is a manifestation of the quaternion identity 

(qm)i(qm) = q{mirh)q. 

We encode the matrix M with the unit quaternion m = Mu — M 2 ii — M 22 j + Mi 2 k, and we will encode and 
Qi with the unit quaternion q = —a + /3i + yj -I- Sk. The Hopf map is visible in the relationship: 

mirh = (j^Mh -\- IW21) ~ (-^12 T -^22)) * T 2(AfiiiWi2 -I- M 21 M 22 ) j + 2 [MiiM 22 ~ M 21 M 12 ) k 
= [i j k] Hopf(M). 

The matrix Q 3 may be found in the computation of q(iX + jY + kZ)q whose {i,j,k) parts are the components 
of Q 3 [X Y Zf. 

The matrix Q 4 may be found in the computation of the quaternion product w = qm by writing W = 

This reverses how m was formed from M. We then create Q 4 from W = —Q 4 M by flattening M and W. 

2.6 The Hemisphere: Positions of special triangles 

We conclude the nonrandom section of this paper with a summary of the correspondence between position on the 
hemisphere and special triangles: 

• Lines of Latitude: Triangles of equal area. 

• Lines of Longitude: Triangles resulting from a fixed rotation of the reference triangle followed by scaling along 
the X and y axes. 

• Circular arc in a plane through an edge of the little equilateral triangle E: Triangles with constant angle. 
(Around the perimeter of the circle, most points correspond to the degenerate triangles with angles (0,0,7r). 
At the vertices of E, we obtain all the triangles {9,tt — 9,0). 

• Vertical circular arcs as above: Right triangles. 

• North Pole: Equilateral triangle. 

• Longitudes at multiples of tt/S: Isosceles triangles 

• The Equator: Singular triangles. 

• Inner Spherical Triangle: Acute triangles 

• The six fold symmetries: The six permutations of unequal a, b, c. 

• Where is the lower hemisphere?: Replace detM with — detM. Every ordered triangle is “double covered.” A 
signed SVD with rotations U and V would work well. 

3 The Normal Distribution generates random shapes 

Section 2 concentrated on the relationship between triangles, the hemisphere and disk, and 2x2 matrices through 
the SVD. Nothing was random. Here in Section 3, randomness comes into its own with the normal distribution as a 
natural structure. 



21 



3.1 The Normal Distribution 


The normal distribution has a very special place in mathematics. Several decades of research into Random Matrix 
Theory have shown that exact analytic formulas are available for important distributions when the matrix entries 
arise from independent standard Gaussians. Non-Gaussian entries rarely enjoy the same beautiful properties for 
finite matrices, though recent “universality” theorems show convergence for infinite random matrix theory. 

We might anticipate, if only by analogy, that random triangles arising from independent Gaussians might also 
enjoy special analytic properties compared to other distributions. Section 3.2 reveals that random triangles generated 
from a Gaussian random matrix are uniformly distributed on the hemisphere. 

Four Gaussian degrees of freedom from six: 

Figure 2 illustrates 1000 random triangle shapes. As we saw in Section 1.3, six independent standard Gaussians 
may be used to describe the 2x3 array of vertex coordinates with six degrees of freedom. 

Nonetheless, if we translate the centroid to the origin, or a vertex to the origin, or use the 2x3 edge matrix E, 
or the matrix from the edge viewpoint or the matrix My from the vertex viewpoint, only four degrees of freedom 
remain. 

Consider the 2x3 matrices T, the origin centered vertex matrix, or E the edge matrix. Up to scaling, both of 
these matrices consist of elements that are standard normals conditioned on a zero column sum. This may be argued 
or computed. The argument is that Gaussians are linear, and the symmetry of the situation requires that no one 
vertex or edge can be special. 

One can check readily that centering six standard normals at the origin produces two independent rows each with 
the singular covariance matrix 



The matrix A^A is the projection matrix onto the plane x -\- y -\- z = 0. The edge matrix E obtained by taking 
differences of columns can be readily verified to have covariance matrix 3A^A. As shape information removes scale 
information, this is irrelevant to the paper, and we often tend to think of T and E as providing dual views of the 
same construction. 

If one forms the matrix My = TA^, the covariance matrix for each row of two elements is now A(A^A)A^ = 12 - 
We then obtain the nice conclusion that My consists of four independent standard normals. Similarly, My before 
scaling would thus consist of four independent normals with variance 3. Once again we remark that scaling is 
irrelevant to shape theory results justifying our normalization of choosing M to have sum of squares 1. 


3.2 Uniform distribution on the hemisphere 

Proposition 4. These equivalent representations describe the uniform measure on the hemisphere with radius 1/2: 


Longitude 

(j) uniform on [0, 27r) is independent from: 

Height 

uniform on [0,1/2] 

Latitude 

density cos(A)dA on [0,7r/2] 

r 

density 4r(l — on [0,1/2] 

det M\ 

uniform on [0,1/2] 

Triangle area 

uniform on [0, l/-\/48] 
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Proof: The latitude and longitude joint density cos(A)dAd(/) is the familiar volume element on the sphere, where 
latitude A is measured from the equator rather than the zenith. Archimedes of Syracuse (born -287, died -212) 
knew the height formula through his hat box theorem: on a sphere, the surface area of a zone is proportional to its 
height. The area comes from Hero (c. 10-70) or Heron of Alexandria’s formula, believed already known to Archimedes 
centuries earlier. 

Lemma 5. Exponential distributions The sum of squares of two independent standard normals (X 2 ) *5 ex¬ 
ponentially distributed with density \e~^^^. If 61,62 are independent random variables with identical exponential 
distribution then eij{ei 62 ) is uniformly distributed on [ 0 , 1 ]. 

Proof: These facts are well known. The generalization to n exponentials is a popular way to generate uniform 
samples Xi = Cij^ti on the simplex Xi = l,Xi>Q (cf. [Ill Proposition 1 ]). 


We now turn to the beautiful result known to David Kendall and his collaborators: 


Theorem 6. Triangles generated from six independent normals ( xy coordinates of three vertices or edges) correspond 
to points distributed uniformly on the hemisphere. 

The corresponding M is a 2 x 2 matrix of independent standard normals. 

Proof 1 (Exponentials): We first consider the height: 


detM/|!M||| = 


ad — be 
-\- c^ df 


(( 


l-( 

[ a—d 


/ b-c 

\V2, 

)“) 

2 

((w)' + (7f) 

2 

+ 1 

f a—d \ 

KV 2 ) 


b-c\ 

.V 2 ) 

0 


The numerator is the difference between two exponentially distributed random variables (being the sums of 
squares of two independent standard normals.) Then by Lemma 5, 


61 - 62 _ 61 1 

2(61 -f 62) 61 -f 62 2 

is uniform on [—1/2,1/2]. The absolute value is then uniform on [0,1/2] as desired. 

If we now turn to the “longitude,” the right singular vectors provide the uniform distribution owing to the right 
orthogonal invariance of M. 

Proof 2 (Random Matrix Theory and condition numbers): 

The singular values (Ti,cr 2 and singular vectors V for M=randn(2,2) are independent. H is a rotation with angle 
uniformly distributed on [ 0 , 7 r). 

The distribution of the condition number k is important. It may be found in [31 Eq. 2.1] or jU Eq. 14]. Restated, 
the probability density of n for a random 2x2 matrix of iid normals is . Equivalently, the distribution 

of is uniform on [ 0 , 1 / 2 ]. 

Proof 3 (Hopf Map) 

We know that M is uniformly distributed on the sphere 1 = -|- -|- M /2 + MI 2 ■ We seek to explain why 

Hopf(M) is uniform on the sphere 1 = x'^ -h y^ -h . The mathematician’s favorite proof would have the latter 
inherited from the former. We want uniformity of Hopf(M) to be a shadow of uniformity of M, but the Hopf map 
is not a linear projection. 
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For any 3x3 rotation matrix Q 3 , we constructed in Section 2.5.6 a 4x4 rotation matrix Q 4 , such that 

Hopf((54M) = (53Hopf(M), 

for any fixed matrix M. If M is random and uniformly distributed on 1 = + M 21 + M 12 + M 22 , then of course 

Q 4 M is too. 

What does this say about the distribution of Hopf(M)? For any choice of Q 3 , the distribution of QsRopflM) is 
the distribution of Hopf(Q 4 M), which is the same as the distribution of Hopf(M). The distribution of Hopf(M) is 
invariant under any fixed 3x3 rotation. It must be the uniform distribution on + y'^ + = 1. 

We encourage the reader to follow this - it is truly a beautiful argument. 

3.3 What is the probability that a random triangle is acnte or obtnse ? 

Put Lewis Carroll’s question in the context of the normal distribution. We then obtain the result known to Portnoy 
and Kendall and collaborators: 

Theorem 7. A random triangle from the uniform distribution in shape space has squared side lengths uniformly 
distributed on [0,2/3]. The probability is 1/4 that this triangle is acute. 

Proof: The edges of the triangle are the lengths of the three columns of MA. Taking the third column for 
convenience, the squared side length is A§ = 2/3 times + MI 2 ). By Lemma 5, this is the uniform distribution 
on [0,2/3]. 

Edges that satisfy = I give a right triangle when = 1/2. The triangle is obtuse when > 1/2, 

which makes c the longest side. The probability that a particular angle is obtuse is then (2/3 — l/2)/(2/3) = 1/4. 
The probability that any angle is obtuse is then 3/4 (at most one can be obtuse!). Then 1/4 is the probability that 
all are acute. □ 

3.4 Triangles in n dimensions 

There is an obvious generalization of Figure 5 to higher dimensions. Let M he & random matrix of independent 
standard norms with n rows and 2 columns. Then MA is a random triangle shape in n dimensions. We can think 
of this as random vertices centered at the origin, or random edges that close to a proper triangle. There is a further 
generalization 

The same argument as for the plane shows a squared side length has the distribution (2/3)Beta(n/2, n/2). In M" 
the probability of an obtuse triangle is3(l — /(|,|,§)) where I denotes the incomplete beta function. This can be 
evaluated in MATLAB by 

3*(betaine(3/4,n/2,n/2,’upper’)) 

or in Mathematica by 

3*(l-BetaRegularized[3/4, n/2, n/2]). 

This probability is also computed by Eisenberg and Sullivan [5]. They note that for larger n, it is increasingly 
likely that a random triangle is acute. The probability of an obtuse triangle is 10% in and 1% in It is 
nearly 0 . 1 % in 

For arbitrary n the distribution is no longer uniform on the hemisphere. 

Connecting multivariate statistics and numerical analysis, the same result is innocently hidden as an exercise 
in Wishart matrix theory. It is an unlikely connection. This “square root ellipticity statistic” may be found in 
Exercise 8.7(b) of [131 P-379]. It states that P < As n increases, this reduces the probability 

near the equator and adds weight near the poles. Acute triangles become more probable. 
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3.5 An experimental investigation that revealed the distribntion’s “shadow” 


Early in our own investigation, we plotted the normalized squares of randn(2, 2) * A as barycentric coordinates. We 
found very quickly that the probability of an acute triangle is 1/4, and triangles naturally fill a hemisphere with a 
uniform distribution. Here is a MATLAB code and a picture that tells so much in so little space that we could not 
resist sharing. A few interesting things to note: Line 10 projects the three normalized squared edges to using 
A"”". Line 9 compares to to decide acute or obtuse. (This is more efficient than computing angles.) Ten 

thousand trials gave 25.37% acute triangles. Ten million trials gave 24.99% acute triangles. 

Readers may wish to recreate this road to discovery of the uniform distribution on the hemisphere: Histogram 
the radii of the points, guess the measure, and then realize the uniformity. We guessed that uniformity by fitting 
the density f{x) = 4a:/Vl — = —1^\/1 — 4a;^ which is the shadow of the uniform distribution on the hemisphere 

(differential form of Archimedes’ hat box theorem). Then came proofs. 


1 t=20000; % trials 

2 v=zeros(3,t); % computed squared edges 

3 D=diag(sqrt([2 6]))\[1 -1 0;1 1 -2]; 

4 for i=l:t 

5 V ( : , i) =sum ( (randn (2) *D) . ^^2 ) ; 

6 V( : , i) =v(:,i)/sum (v ( : , i) ) ; 

7 end 

8 figure(1);hold off 

9 obtuse=any(v>0.5) ; 

10 twod=D*v*sqrt(3/2) ; 

11 plot(twod(1,obtuse),twod(2, obtuse) , ' . k ' ) 



Figure 10: Random a^,6^,c^ are computed on line 5 and normalized on line 6. The obtuse triangles are identified 
on line 9. Line 10 projects the plane x + y + z = 1 onto Line 11 plots every obtuse triangle as a point in this 
plane. The acute triangles form the inner triangle as in Figure 3 and the obtuse triangles complete a disk of radius 
1/2. The hemisphere is viewed from above. 



Figure 11: Line 12 computes the sample probability of acute triangles. Lines 13 through 18 histogram all the radii 
against the guess (curved line) that the points are uniformly distributed on the hemisphere. 
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3.6 Uniform shapes versus uniform angles 

Two uniform distributions, on the hemisphere and on angle space A + B -\- C = 180°, gave the same fraction | of 
obtuse triangles. We are not aware of a satisfying theoretical link. Portnoy m Section 3] philosophizes about “the 
fact that the answer 3/4 arises so often”. To emphasize the difference between these distributions, we report on a 
numerical experiment and a theoretical density computation to understand the angular distribution. 

Our first figure might be called lOQ 000 triangles in 100 bins. The three angles divided by tt are barycentric 
coordinates in the figure. With four bins the triangles would appear uniform. With 100 bins we see that they are 
anything but. A uniform distribution would have 1000 triangles per bin. 

100,000 triangles in 100 bins 



C u:- .VySt o 

Green) 10-25% Percentile 
Blue) 25 - 75 % Percentile 
Magenta) 75-90% Percetile 
Top 10% 


Figure 12: 100, 000 triangles in 100 bins: Triangles selected uniformly in shape space are not uniform in angle space. 

The underlying theoretical distribution involves the Jacobian from angle space to the “squared edges” disk. Then 
one must apply a second Jacobian to the hemisphere, and finally invert. Since the second part is standard, we will 
sketch the first part of the argument. 
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Suppose a + ^ + 7 = 1. Considering the Law of Sines, define 


s = 



^ I sin^ na 
— ( sin^ Tr/3 

sin^ 7r7 


with a = sin^ na + sin^ tt/? + sin^ ttj. With some calculus we can show that 


/ da \ 

ds = TT J I dp I , 

V dl ) 


T ( sin27ra \ 

with J = (—— ^-) and p = sin27r/? . The Jacobian determinant that transforms the angle space 

y sin 27r7 j 

to squared edge space is proportional to det (AJA"'"). The remaining step, omitted here, is the Jacobian to the 
hemisphere. As the experiment begins with triangles uniform on the hemisphere, it is the inverse Jacobian that we 
use in our plot. 

Figure 12 is a Monte-Carlo experiment. Figure 13, shows the theoretical Jacobian. 



Figure 13: Angle density reveals the nonuniformity of angles, when triangles are picked uniformly from shape space. 
The least likely angles are in black, cyan and yellow with p < 0.1, p < 0.25, and p < .4. The most likely angles are 
in blue, magenta, and red with p > 0.6, p > 0.75, and p > 0.90. Green is .4 < p < .6. 
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4 Further Applications 

A corpus of pictures becomes a set of shapes, and we study a particular feature. A random choice then becomes a 
random matrix and the entirety of random matrix theory can be considered for shape applications. 

We believe that the timing is right to realize the dream of shape theory. In our favor are 

• New theorems in random matrix theory that can immediately apply to shape theory 

• Modern computational power making shapes accessible (and Monte Carlo simulations) 

• The technology of multivariate statistical theory. We can compute hypergeometric functions of matrix argument 
and zonal polynomials m- Muirhead m was 30 years ahead of his time. 

There is much room for research in this area. In two dimensions, the exact condition number distribution for the 
sphericity test confirms that the shapes come from the standard normal distribution. 

It is very possible that uniformity may not be a great measure for real problems on shape space. Similar issues 
arise for random matrices and random graphs. It is a mistake to think of a random object as being “any old object.” 
Usually random objects have special properties of their own, much like a special matrix, or graph, or shape. 

Here are random shapes and convex hulls. The general technique is multiplying a matrix of standard normals by 
a A„ from the Helmert matrix. It is a bit further afield (but interesting) to ask for the average number of edges of 
the convex hull. For literature in this direction see [ini Chapter 8]. 



k = 100 


# 4 ^ A < 9 - » 


Figure 14: Ten random shapes in 2d taken from the uniform shape distribution with fc = 3, 5, 20,100 points. 
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Figure 15: Random tetrahedra in 3d (m = 3 and fc = 4) 
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Figure 16: Random Gems: Convex hulls of 100 random points in 3d (m = 3 and k = 100) 


4.1 Tests for Uniformity of Shape Space 

An object formed from k points in K™ may be encoded in an m x k matrix X. The shape is therefore encoded in 
an TO X (fc — 1) preshape matrix Z = /||AA^||^, where Aj. is the Helmert matrix in equation Q. The uniform 

distribution on shapes may be thought of as the distribution obtained when 


X = randn(m,k-l) * A^, 


so that X is the product of an to x (k— 1) matrix of iid standard Gaussians and A^. In the preshape, Z is normalized 
by its own Frobenius norm, 

Z — randn(m,k-l)/|| • H^;’. 

Figure 15 plots random tetrahedra in this way. 

In [2], Chikuse and Jupp propose a statistical test on Z for uniformity. From samples Z\^... ^Zt they calculate 


S = 


(k — l){m{k — 1) + 2) 


t trace({iX;Li^'^^- ) 


As t —>■ oo they approximate S with x\h-i)(k+2)/2- Corrections are proposed for finite t as well. 

Other tests are easy to construct. For example, much is known about the smallest singular value of the random 
matrix Z. When Z is square and t > m, the density of l/cri„in(-^) is exactly 


- 1 ) 






2F1 




vn? + TO 
2 
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Non-square cases can also be handled by a combination of known techniques. 

With the exact density, one can compute the smallest singular value of the samples and perform goodness-of-fit 
tests (such as Kolmogorov-Smirnov). Broadly speaking one might choose to test for the orthogonal invariance of the 
singular vectors. Since Z * Xm(fc-i)is normally distributed, tests based on Wishart matrices are natural. 

4.2 Northern Hemisphere Map 

Perhaps because we had the technology, we mapped the northern hemisphere (shape theory hemisphere) to angle 
space (Figure 17.) Barycentric coordinates correspond to the angles divided by tt. The resulting picture appears in 
the figure that follows. The middle triangle consists of the “acute” points. 



Acute Territon 

\ 


Figure 17: Northern hemisphere as triangles in shape space, mapped to “angle” space. 

We would like to thank Wilfrid Kendall, Mike Todd, and Eric Kostlan for their insights. The first author 
acknowledges NSF support under DMS 1035400 and DMS 1016125. The second author acknowledges NSF support 
under EFRI 1023152. 
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